Deep-sequencing identification of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

Multilocus sequence analysis and rpoB sequencing of Mycobacterium abscessus (sensu lato) strains

Mycobacterium abscessus, Mycobacterium bolletii, and Mycobacterium massiliense (Mycobacterium abscessus sensu lato) are closely related species that currently are identified by the sequencing of the rpoB gene. However, recent studies show that rpoB sequencing alone is insufficient to discriminate between these species, and some authors have questioned their current taxonomic classification. We studied here a large collection of M. abscessus (sensu lato) strains by partial rpoB sequencing (752 bp) and multilocus sequence analysis (MLSA). The final MLSA scheme developed was based on the partial sequences of eight housekeeping genes: argH, cya, glpK, gnd, murC, pgm, pta, and purH. The strains studied included the three type strains (M. abscessus CIP 104536(T), M. massiliense CIP 108297(T), and M. bolletii CIP 108541(T)) and 120 isolates recovered between 1997 and 2007 in France, Germany, Switzerland, and Brazil. The rpoB phylogenetic tree confirmed the existence of three main clusters, each comprising the type strain of one species. However, divergence values between the M. massiliense and M. bolletii clusters all were below 3% and between the M. abscessus and M. massiliense clusters were from 2.66 to 3.59%. The tree produced using the concatenated MLSA gene sequences (4,071 bp) also showed three main clusters, each comprising the type strain of one species. The M. abscessus cluster had a bootstrap value of 100% and was mostly compact. Bootstrap values for the M. massiliense and M. bolletii branches were much lower (71 and 61%, respectively), with the M. massiliense cluster having a fuzzy aspect. Mean (range) divergence values were 2.17% (1.13 to 2.58%) between the M. abscessus and M. massiliense clusters, 2.37% (1.5 to 2.85%) between the M. abscessus and M. bolletii clusters, and 2.28% (0.86 to 2.68%) between the M. massiliense and M. bolletii clusters. Adding the rpoB sequence to the MLSA-concatenated sequence (total sequence, 4,823 bp) had little effect on the clustering of strains. We found 10/120 (8.3%) isolates for which the concatenated MLSA gene sequence and rpoB sequence were discordant (e.g., M. massiliense MLSA sequence and M. abscessus rpoB sequence), suggesting the intergroup lateral transfers of rpoB. In conclusion, our study strongly supports the recent proposal that M. abscessus, M. massiliense, and M. bolletii should constitute a single species. Our findings also indicate that there has been a horizontal transfer of rpoB sequences between these subgroups, precluding the use of rpoB sequencing alone for the accurate identification of the two proposed M. abscessus subspecies.

Ambiguous nucleotide calls from population-based sequencing of HIV-1 are a marker for viral diversity and the age of infection

The time passed since the infection of a human immunodeficiency virus (HIV)-infected individual (the age of infection) is an important but often only poorly known quantity. We assessed whether the fraction of ambiguous nucleotides obtained from bulk sequencing as done for genotypic resistance testing can serve as a proxy of this parameter.

Targeted next generation sequencing as a diagnostic tool in epileptic disorders

Epilepsies have a highly heterogeneous background with a strong genetic contribution. The variety of unspecific and overlapping syndromic and nonsyndromic phenotypes often hampers a clear clinical diagnosis and prevents straightforward genetic testing. Knowing the genetic basis of a patient's epilepsy can be valuable not only for diagnosis but also for guiding treatment and estimating recurrence risks.

Identification of the bovine Arachnomelia mutation by massively parallel sequencing implicates sulfite oxidase (SUOX) in bone development

Arachnomelia is a monogenic recessive defect of skeletal development in cattle. The causative mutation was previously mapped to a approximately 7 Mb interval on chromosome 5. Here we show that array-based sequence capture and massively parallel sequencing technology, combined with the typical family structure in livestock populations, facilitates the identification of the causative mutation. We re-sequenced the entire critical interval in a healthy partially inbred cow carrying one copy of the critical chromosome segment in its ancestral state and one copy of the same segment with the arachnomelia mutation, and we detected a single heterozygous position. The genetic makeup of several partially inbred cattle provides extremely strong support for the causality of this mutation. The mutation represents a single base insertion leading to a premature stop codon in the coding sequence of the SUOX gene and is perfectly associated with the arachnomelia phenotype. Our findings suggest an important role for sulfite oxidase in bone development.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Real-time quantitative broad-range PCR assay for detection of the 16S rRNA gene followed by sequencing for species identification

Here we determined the analytical sensitivities of broad-range real-time PCR-based assays employing one of three different genomic DNA extraction protocols in combination with one of three different primer pairs targeting the 16S rRNA gene to detect a panel of 22 bacterial species. DNA extraction protocol III, using lysozyme, lysostaphin, and proteinase K, followed by PCR with the primer pair Bak11W/Bak2, giving amplicons of 796 bp in length, showed the best overall sensitivity, detecting DNA of 82% of the strains investigated at concentrations of < or =10(2) CFU in water per reaction. DNA extraction protocols I and II, using less enzyme treatment, combined with other primer pairs giving shorter amplicons of 466 bp and 342 or 346 bp, respectively, were slightly more sensitive for the detection of gram-negative but less sensitive for the detection of gram-positive bacteria. The obstacle of detecting background DNA in blood samples spiked with bacteria was circumvented by introducing a broad-range hybridization probe, and this preserved the minimal detection limits observed in samples devoid of blood. Finally, sequencing of the amplicons generated using the primer pair Bak11W/Bak2 allowed species identification of the detected bacterial DNA. Thus, broad-spectrum PCR targeting the 16S rRNA gene in the quantitative real-time format can achieve an analytical sensitivity of 1 to 10 CFU per reaction in water, avoid detection of background DNA with the introduction of a broad-range probe, and generate amplicons that allow species identification of the detected bacterial DNA by sequencing. These prerequisites are important for its application to blood-containing patient samples.

Evaluation of the colorimetric VITEK 2 card for identification of gram-negative nonfermentative rods: comparison to 16S rRNA gene sequencing

Ninety strains of a collection of well-identified clinical isolates of gram-negative nonfermentative rods collected over a period of 5 years were evaluated using the new colorimetric VITEK 2 card. The VITEK 2 colorimetric system identified 53 (59%) of the isolates to the species level and 9 (10%) to the genus level; 28 (31%) isolates were misidentified. An algorithm combining the colorimetric VITEK 2 card and 16S rRNA gene sequencing for adequate identification of gram-negative nonfermentative rods was developed. According to this algorithm, any identification by the colorimetric VITEK 2 card other than Achromobacter xylosoxidans, Acinetobacter sp., Burkholderia cepacia complex, Pseudomonas aeruginosa, and Stenotrophomonas maltophilia should be subjected to 16S rRNA gene sequencing when accurate identification of nonfermentative rods is of concern.

Mutation screening of the medium-chain acyl-CoA dehydrogenase (MCAD) and the ornithine transcarbamylase (OTC) genes by multiplex PCR amplification and sequencing

BACKGROUND: Sequencing based mutation screening assays of genes encompassing large numbers of exons could be substantially optimized by multiplex PCR, which enables simultaneous amplification of many targets in one reaction. In the present study, a multiplex PCR protocol originally developed for fragment analysis was evaluated for sequencing based mutation screening of the ornithine transcarbamylase (OTC) and the medium-chain acyl-CoA dehydrogenase (MCAD) genes. METHODS: Single exon and multiplex PCR protocols were applied to generate PCR templates for subsequent DNA sequencing of all exons of the OTC and the MCAD genes. For each PCR protocol and using the same DNA samples, 66 OTC and 98 MCAD sequence reads were generated. The sequences derived from the two different PCR methods were compared at the level of individual signal-to-noise ratios of the four bases and the proportion of high-quality base-signals. RESULTS: The single exon and the multiplex PCR protocol gave qualitatively comparable results for the two genes. CONCLUSIONS: Many existing sequencing based mutation analysis protocols may be easily optimized with the proposed method, since the multiplex PCR protocol was successfully applied without any re-design of the PCR primers and other optimization steps for generating sequencing templates for the OTC and MCAD genes, respectively.

The Mycoplasma conjunctivae genome sequencing, annotation and analysis

BACKGROUND: The mollicute Mycoplasma conjunctivae is the etiological agent leading to infectious keratoconjunctivitis (IKC) in domestic sheep and wild caprinae. Although this pathogen is relatively benign for domestic animals treated by antibiotics, it can lead wild animals to blindness and death. This is a major cause of death in the protected species in the Alps (e.g., Capra ibex, Rupicapra rupicapra). METHODS: The genome was sequenced using a combined technique of GS-FLX (454) and Sanger sequencing, and annotated by an automatic pipeline that we designed using several tools interconnected via PERL scripts. The resulting annotations are stored in a MySQL database. RESULTS: The annotated sequence is deposited in the EMBL database (FM864216) and uploaded into the mollicutes database MolliGen http://cbi.labri.fr/outils/molligen/ allowing for comparative genomics. CONCLUSION: We show that our automatic pipeline allows for annotating a complete mycoplasma genome and present several examples of analysis in search for biological targets (e.g., pathogenic proteins).

Genotyping of Francisella tularensis strains by pulsed-field gel electrophoresis, amplified fragment length polymorphism fingerprinting, and 16S rRNA gene sequencing

We evaluated three molecular methods for identification of Francisella strains: pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphism (AFLP) analysis, and 16S rRNA gene sequencing. The analysis was performed with 54 Francisella tularensis subsp. holarctica, 5 F. tularensis subsp. tularensis, 2 F. tularensis subsp. novicida, and 1 F. philomiragia strains. On the basis of the combination of results obtained by PFGE with the restriction enzymes XhoI and BamHI, PFGE revealed seven pulsotypes, which allowed us to discriminate the strains to the subspecies level and which even allowed us to discriminate among some isolates of F. tularensis subsp. holarctica. The AFLP analysis technique produced some degree of discrimination among F. tularensis subsp. holarctica strains (one primary cluster with three major subclusters and minor variations within subclusters) when EcoRI-C and MseI-A, EcoRI-T and MseI-T, EcoRI-A and MseI-C, and EcoRI-0 and MseI-CA were used as primers. The degree of similarity among the strains was about 94%. The percent similarities of the AFLP profiles of this subspecies compared to those of F. tularensis subsp. tularensis, F. tularensis subsp. novicida, and F. philomiragia were less than 90%, about 72%, and less than 24%, respectively, thus permitting easy differentiation of this subspecies. 16S rRNA gene sequencing revealed 100% similarity for all F. tularensis subsp. holarctica isolates compared in this study. These results suggest that although limited genetic heterogeneity among F. tularensis subsp. holarctica isolates was observed, PFGE and AFLP analysis appear to be promising tools for the diagnosis of infections caused by different subspecies of F. tularensis and suitable techniques for the differentiation of individual strains.